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Could AGN Outbursts Transform Cool Core Clusters? 
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ABSTRACT 

The origin of the bimodality in cluster core entropy is still unknown. At the same time, 
recent work has shown that thermal conduction in clusters is likely a time-variable 
phenomenon. We consider if time-variable conduction and AGN outbursts could be 
responsible for the cool-core (CC), non cool-core (NCC) dichotomy. We show that 
strong AGN heating can bring a CC cluster to a NCC state, which can be stably 
maintained by conductive heating from the cluster outskirts. On the other hand, if 
conduction is shut off by the heat-flux driven buoyancy instability, then the cluster 
will cool to the CC state again, where it is stabilized by low-level AGN heating. Thus, 
the cluster cycles between CC and NCC states. In contrast with massive clusters, we 
predict the CC/NCC bimodality should vanish in groups, due to the lesser role of 
conductive heating there. We find tentative support from the distribution of central 
entropy in groups, though firm conclusions require a larger sample carefully controlled 
for selection effects. 
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1 INTRODUCTION 

There is a striking observational bimodality in the prop- 
erties of galaxy cluster cores, which can be broadly sep- 
arated into two types: cool-core (CC) and non cool-core 
(NCC) clusters. The former are defined to have temper- 
ature profiles which decline significantly toward the cen- 
ter. Measuring the relative abundance of the two types is 
somewhat hampered by selection effects (CC clusters have 
strongly peaked X-ray emission profiles which are more eas- 
ily detected), but surveys indicate that roughly half of all 
clusters are in each category (e.g., I Chen et al. [2007). The 
logarithmic slope of the entropy profile is bimodal, with 
CC/N CC clusters having ste eper /shallower slopes respec- 
tively (jSanderson et all [2009). The distribution of core en- 
tropy also appears to be bimodal in clusters, with population 
peaks at K ~ 15keVcm 2 and K ~ 150 ke V cm 2 and a dis- 
tinct gap between K ~ 30 — 50 keV cm 2 (|Cavagnolo et all 
2009). This in turns feeds into star formation and AGN 
properties: Ha and radio emission from the central brightest 
cluster galaxy are much more pronounced when the cluster's 
core falls below an ent ropy threshold of K < 30keVcm -2 
l|Cavagnolo et al.ll2008l ), and a major ity (~ 70%) of CC clus- 
ters host radio sources (lBurnslll990l ). 

Clearly, unravelling the origin of this dichotomy-for 
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which there is no widely accepted explanation-could poten- 
tially yield great insight into cluster thermodynamics. Merg- 
ers have been considered a prime candidate for transform- 
ing CC to NCC systems, given the frequency of mergers in 
a hierarchical CDM cosmology, as well as the large amount 
of energy (as much as ~ 10 64 erg) in mergers, well in ex- 
cess of that required. However, detailed simu lations have 
not b orne this expectation out. For instance, IPoole et alj 
(2008) find that CC systems are remarkably robust and only 
disrupted in direct head-on or multiple collisions; even so, 
the resultin g warm core state is only transient. To date, 
I Burns et al.l l|2008l ) present the only set of simulations where 
NCC clusters are produced via mergers. For this to hap- 
pen, nascent clusters must experience major mergers early 
on which destroyed embryonic CCs and prevented their ref- 
ormation. Note that these simulations do not incorporate 
mechanisms (such as AGN feedback) to stop a cooling catas- 
trophe; furthermore, the relatively low numerical resolu- 
tion (15.6 h~ kpc) may preclude firm concl usions about core 
struct ure and evolution. As an alternative, iMcCarthv et alj 
(2008) suggested that early pre- heating prior to cluster col- 
lapse could explain the lack of low entropy gas in NCC sys- 
tems, which receive higher levels of preheating compared to 
CC systems. A possible concern in such scenarios is whether 
one can pre-heat the ICM to a high adiabat and yet re- 
tain sufficient low entropy gas i n lower mass halo s to ob- 
tain a realistic galaxy population dBower et al 1l2008l ; sec also 
lOh fc Benson! l2003l : IScannapieco fc Ohll2004 ). More impor- 
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tantly, many N CC clusters also have a short central cooling 
time (~ 1 Gvr: ISanderson et alJ 12006). and it is not clear 
why radiative cooling should not erase memory of the ini- 
tial preheating episode. Both of these hypotheses focus on 
'nature', or initial conditions, in the form of an early major 
merger or preheating, in determining whether a cluster is 
CC or NCC. This hints at a fine-tuning problem: why are 
initial conditions such that roughly equal numbers of CC 
and NCC systems appear? In addition, there appears to be 
substantial differences in the metallicity profiles of NCC and 
CC systems, at least in groups (see discussion). Simulations 
show that metall icity profiles are r emarkably stable to sub- 
sequent mergers ( Poole et al.ll2008h . 

Recently, in iGuo et aUl|200a , hereafter GOR08), we 
conducted a global Lagrangian stability analysis of clusters 
in which cooling is balanced by AGN heating and thermal 
conduction. This offered an alternative promising explana- 
tion, based on 'nurture', or physical processes occurring in 
the ICM in its recent past. Our analysis showed that glob- 
ally stable clusters could only exist in two forms: (1) cool 
cores stabilized by both AGN feedback and conduction, or 
(2) non-cool cores stabilized primarily by conduction^. In- 
termediate temperature profiles typically lead to globally 
unstable solutions, which would then quickly evolve to ei- 
ther CC or NCC states. In GOR08, we speculated that 
these two categories of clusters might even represent dif- 
ferent stages of the same object. The importance of thermal 
conduction on global scales obviously depends on the large 
scale structure of the cluster magnetic fields. Recent cal- 
culations suggest that thermal conduction of heat into the 
cluster core can be self-limiting: in cases where the tempera- 
ture decreases in the direction of gravity, a buoyancy insta- 
bility (the heat flux driven buoyancy instability, hereafter 
HBI) sets in which re-orients a radial magnetic field to be 
larg ely transverse, shutting off conduction to the cluster cen- 
ter jQuataert 2008; iParrish fe Quataerdlioosl : iParrish" et al.l 
l2009l : lBogdanovic et alj|2009h. In GOR08 , we speculated (as 
subsequently did lBogdanovic et al.ll2009l ) that powerful out- 
bursts from a central AGN might counteract the HBI by 
disturbing the azimuthal nature of the magnetic field, thus 
enabling thermal conduction. In particular, the following 
scenario could arise: as conductivity falls, gas cooling and 
mass inflow will increase, triggering AGN activity. The ris- 
ing buoyant bubbles may re-orient the magnetic field to be 
largely radial again, increasing thermal conduction and re- 
ducing mass inflow, shutting off the AGN until the HBI sets 
in once again. If AGN heating and/or thermal conduction 
during their 'on' states are strong enough to heat the CC 
cluster to the NCC state, the cluster could then continuously 
cycle between cool-core (AGN heating dominated) and non 
cool-core (conduction dominated) states. 

The goal of this paper is to perform a very simple fea- 
sibility study for such a scenario, motivating future, more 
detailed work. While the transformation of NCC to CC clus- 
ters via radiative cooling can be easily accomplished (as 



1 A recent study of a Chandra cluster sample has similarly found 
that while thermal conduction appears to be sufficient to stabilize 
NCC clusters, CC clusters appear to form a di stinct population 
in which additional feedback heating is required ( Sanderson ct al. 
|200S|) 



seen, for instance, in Fig. 13 of IParrish et all ((2009)), the 
possible transformation of CC to NCC systems via AGN 
outbursts or time variable conduction has not been demon- 
strated. It has been conjectured before that strong AGN 
outbursts (the most extreme examples of which are Hydra 
A and MS 0435+7241) could permanently modify co re en- 
tropy (|Kaiser fe Binnevl 120031 : IVoit fe Donahue! I2005T ). But 
there have been no explicit calculations, and indeed sugges- 
tions that a CC t o NCC transformation would be energet- 
ically prohibitive l|McCarthv et al]|2008h . In this paper, we 
perform explicit calculations to investigate if strong AGN 
outbursts could transform a CC cluster to the NCC state 
( ij3.2jl . Even if energetically subdominant (as they likely are 
for the most massive clusters), AGN could catalyze a dom- 
inant contribution from thermal conduction either due to: 
(i) the strong temperature dependence of thermal conduc- 
tive flux, F tx T 5 / 2 ; (ii) altering magnetic topology as dis- 
cussed above. We study if such effects can indeed allow a 
CC to NCC transformation. It is important to note that 
the ability of rising bubbles or other bulk gas motions to 
globally restructure field geometry and hence thermal con- 
ductivity has not been demonstrated. However, there is im- 
portant circumstantial support from numerical simulations 
of magnetic draping, which show that magnetic fields are 
amplified and m ore ordered in the wak e of moving sub- 
halos or bubbles dRuszkowski et al.l 2007 : Asai et alj 120071 ; 



iDursi fe Pfrommerll2008l ; io'Neill et al.ll2009l ). To isolate the 
relative contribution of AGN outbursts and conduction, we 
also consider models of galaxy groups. The strong tempera- 
ture dependence of thermal conduction implies that conduc- 
tion should be irrelevant in groups, regardless of magnetic 
field geometry. If our explanation for the bimodality in clus- 
ter cores is correct, then such bimodality should disappear 
in galaxy groups. 

We describe our methods in §2, our results in §3, and 
discuss implications in §4. The cosmological parameters used 
throughout this paper are: fl m = 0.3, Qa ~ 0.7, h — 0.7. 
We have rescaled observational results if the original paper 
used a different cosmology. 



2 BASIC ASSUMPTIONS AND SETUP 

We solve the time-dependent hydro dynamic equations usin 
the ZEUS-3D hydrodynamic code (|Stone fe Normanlll99~ 
in its one-dimensional mode; in particular, we have in- 
corporated into ZEUS a background gravitational poten- 
tial, radiative cooling, thermal conduction, conv ection, and 
AGN heating. We used the code similarly in IGuo fe Ohl 
(2008), albeit modified for the case of cosmic-ray heating, 
and gratefully acknowledge Mateusz Ruszkowski for sup- 
plying us with the modified ve rsion of ZEUS described in 
iRuszkowski fc B cgclman ( 2002;, hereafter RB02) , which was 
used as our base code. The model of AGN heatin g we adopt 
here is that of 'effervescent heating' proposed by ([Begelmanl 
l200lt) . simulated in RB02, and also used in the global sta- 
bility model of GOR08. We refer the reader to these papers 
for details of the models. Here we simply summarize several 
modifications and reiterate some important points to note. 

In all of the preceding papers, it was assumed that the 
AGN kinetic luminosity was directly related to the instan- 
taneous mass accretion rate, L = eMc 2 . This in itself was 
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a simplification, given that AGN activity is likely to be in- 
termittent, and not necessarily directly related to the mass 
inflow rate. It was justified on the grounds that AGN in- 
termittency timescales are likely shorter than bubble rise 
times and gas cooling times, and hence that AGN feedback 
can be incorporated in a time-averaged sense. Here, we drop 
this assumption, and instead make the alternative (and per- 
haps more realistic) prescription of directly incorporating 
AGN outbursts and intermittency in the simulations. In our 
model, we assume that AGN is triggered when the central 
gas entropy drops below a critical value (S5 < S cr it, where 
55 is the gas entropy at radius of r = 5 kpc from the cluster 
center). We have also considered other AGN trigger crite- 
ria: for example, the AGN is triggered when the central mass 
inflow rate is larger than a critical value or the central gas 
cooling time is less than a critical value. We found that our 
results are fairly robust to the specific criterion we adopted. 
See Table Q] for the specific AGN trigger criterion adopted 
for each run. 

Once an AGN is triggered, we assume that the AGN 
heating lasts for a duration of order the bubble rise time, 
which is typically comparable to the sound crossing time 
t sc ~ 10 8 noocJ 1000 yr for a radius r ~ lOOrioo kpc and 
sound speed c B ~ 1000c Sj iooo kms -1 (e.g., see table 3 in 
iBirzan et all |2004 )'). Thus we assume that each AGN heat- 
ing episode starts once the AGN trigger criterion is satisfied, 
and lasts for i agn ~ 1 x 10 8 yr. During the outburst, we as- 
sume a energy E agn ~ 10 60 — 10 61,5 erg is liberated. Estimates 
of the work needed to i nflate observed cavi t ies in rich clusters 
yield s E agn ~ 10 60 erg (jBirzan et al.ll2004l : iMcNamara et all 
120051 ), rising by 1-2 orders of magnitude in the most ex- 
treme events such as Hydro A dNulsen et alj|2005l ) and MS 
0735+7421 ()McNamara et al.ll2005r i. By comparison, note 
that a ~ 10 9 Mq black hole which doubles its mass over an 
accretion episode will liberate ~ 10 62 (e/0.05) erg, where e is 
the efficiency in converting rest mass energy to kinetic en- 
ergy. The AGN luminosity during each active AGN heating 
episode is then L = -E agn /t agn ~ 10 44 — 10 46 ergs -1 . After a 
period i agn with active AGN heating, we turn off the AGN 
heating, but continue to monitor the AGN trigger criterion. 
Once it is again satisfied, a new AGN heating episode starts. 
The AGN duty cycle (the period between the triggering of 
two successive AGN episodes) is set by the cooling time 
and is generally larger than t agn . During each outburst, the 
AGN heating rate is determined by the 'effervescent heat- 
ing' model described in RB02, except that here we adopt a 
stronger inner heating cutoff term 1 — e~^ r ^ T °^ , instead of 
1 — e~ r//r ° adopted in RB02, to account for the finite size 
of the central radio source and to avoid overheating at the 
cluster centre. In the rest of this paper, the inner heating 
cutoff radius ro is taken to be 10 kpc. 

In our simulations, the gas entropy usually increases 
with radius, agreeing with observational trends. However, 
when the ICM is heated by a strong AGN outburst, neg- 
ative entropy gradients may appear in some regions for a 
short time period. Thus, convection is also included in our 
calculation; the convective flux -F conv is given by the mixing 
length theory described in RB02. In cluster regions with neg- 
ative entropy gradients, convection is turned on and trans- 
ports thermal gas energy. Similar to RB02, we found that 
convection is not important for the parameters of the mod- 



els presented in this paper: we ran our simulations for the 
same models but without convection and found similar re- 
sults. In the low-density weakly-magnetized plasma (e.g., the 
ICM), anisotropic thermal conduction preferentially along 
the magnetic field lines modifies the usual Schwarzschild 
convectiv e criterion th r ough the magneto-thermal instabil- 
ity fMTI: lBalbus]|200d ; |Parrish fc Stonell2005l , 120071 ). While 
we do not incorpora te magnetic fields or MTI, recent cal- 
culat ions which do l|Bogdanovic et all 120091 : iParrish et al.l 
2009) similarly find it to be unimportant by two orders of 
magnitude. 

We do not conduct self-consis tent MHD s i mula- 
tions of the HBI (as for instance in IParrish et alj 120091 : 
iBogdanovic et aTll2009l ) and so instead employ a simplified 
toy model for the conductivity. Given that the mechanism 
for overcoming the HBI is as yet unknown, we feel that such 
illustrative examples of the possible impact of time- variable 
conductivity is justified. We assume that radial conductivity 
is characterized by the Spitzer conductivity with some time- 
dependent suppression factor /. Initially, the conductivity is 
assumed to be negligible; we then assume that thermal con- 
duction is efficient (see Table [T] for the value of / adopted 
in each run) during each AGN outburst (t agn ~ 1 X 10 8 yr) 
and then is either turned off or decays as ~ e -(*- t o£f)/*HBi ) 
where t s is the end time of the preceding AGN heating 
episode. In practice, we have found that either assumption, 
as well as simulations where the onset of efficient conduc- 
tion lags behind the AGN trigger, all yield very similar re- 
sults with regard to the stability of the CC state. Here the 
HBI growth time £hbi ~ 1.0 x 10 8 yr for the typical cool 
core cluster A1795 and may be much lar ger for non-coo l 
core clusters, since £hbi oc (dlnT/dr) 1 / 2 (|Quataer"t1 120081 : 
IParrish et ai1l200St ). 

Our computational grid extends from r m i n (1 kpc) to 
r max (200 kpc for A1795 and 100 kpc for NGC 4325). In 
order to resolve adequately the inner regions, we adopt 
a logarithmically spaced grid in which (Ar)i+i/(Ar)j = 
(fmax/fmin) 1 , where N is the number of active zones. The 
standard resolution of our simulations presented in this pa- 
per is N = 400; our code has been tested to be numerically 
convergent through simulations with different levels of reso- 
lutions. 

For initial conditions, we assume the ICM to be isother- 
mal at T = T ou t, and solve for hydrostatic equilibrium. We 
assume that at the outer boundary r max , n e (r max ) = n out , 
which is close to the value extrapolated from the observa- 
tional density profile. For boundary conditions, we assume 
that the gas is in contact with a thermal bath of constant 
temperature and pressure at the outer radius, where the 
cooling time exceeds the Hubble time. Thus, we ensure that 
temperature and density of the thermal gas at the outer 
radius are constant. We extrapolate all hydrodynamic vari- 
ables from the active zones to the ghost zones by allowing 
them to vary as a linear function of radius at both the in- 
ner and outer boundaries. The intracluster gas is allowed to 
flow in and out of active zones at both the inner and outer 
boundaries. 
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Figure 1. Time evolution of entropy at r = 5 kpc in various 
models for the cluster Abell 1795 (top) and for the group NGC 
4325 (bottom). A strong AGN outburst and/or a sharp increase 
in conductivity is able to accomplish a transition from a CC to 
NCC state. 

3 RESULTS 

3.1 Stability in the CC state 

Let us begin by considering how clusters can be stabilized 
in the CC state by a combination of AGN and conduc- 
tive heating. We first run simulations for a typical mas- 
sive cluster Abell 1795. The parameters for this cluster 
are M = 6.6 x IP 14 M^ , r a = 460 kpc, r c = r s /20 
( Zakamsk a fc Naravanl 120031 ; also see GOR08 for details) , 
Tout = 6.8 keV, and n out = 0.003 cm -3 . Specific model pa- 
rameters for each run are listed in Table [T] 

Run A1795-1 is a representative simulation. The clus- 
ter is initially in hydrostatic equilibrium with spatially con- 
stant temperature T = T out = 6.8 keV, and then evolves 
via radiative cooling without AGN heating or thermal con- 
duction. The solid line in the top panel of Figure [1] shows 
the time evolution of gas entropy at r — 5 kpc (Ss), which 
decreases gradually in the first 4 Gyr. When the central gas 
entropy drops below 15 keV cm 2 (at t ~ 4 Gyr), an AGN 
outburst with E agn — 2 x 10 60 erg is triggered and lasts for 
iagn = 1 x 10 s yr. During the period of active AGN heating, 
thermal conduction with / = 0.4 is also turned on. As clearly 
shown in Fig. [TJ the cooling catastrophe is quickly averted, 
and the gas entropy increases. After 0.1 Gyr of active heat- 
ing, both AGN heating and thermal conduction are turned 
off and the cluster cools again until the next heating episode 
is triggered. As seen for t ~ 4 — 6 Gyr in Fig. [T] (top panel), 
the ICM entropy executes minor oscillations and the clus- 
ter stays in the cool core state, where radial p rofiles of gas 
temp erature and density fit observational data (|Ettori et al.l 
2002) very well. The radial profiles of entropy oscillating in 
the CC state are shown in Figure [2] (top), where entropy 
profiles are plotted every 0.1 Gyr since t = 4 Gyr. During 



the CC state, the lines are clearly concentrated in the lower 
branch with central gas entropy ~ 10 — 20 keV cm 2 . 

In run A1795-1, the AGN duty cycle is ~ 0.3 Gyr. Dur- 
ing each heating episode, the volume-integrated conductive 
heating energy is around 2£ , agn , i.e., the conductive heating 
is comparable to AGN heating. We have also done similar 
calculations for the cluster A2199 (T ou t = 4.6 keV) and the 
group NGC 4325 (T ou t = 1 keV), and found that conduc- 
tive heating is an order of magnitude smaller than AGN 
heating in the former and becomes negligible in the lat- 
ter. The results for NGC 4325 are also presented in this 
paper for comparison. The parameters for this group are 
M = 1.1 x 10 13 M Q , r s = 78.3 k pc, r c = kpc, T^ t = 1 
keV, and n out = 7.27 x 10~ 4 cm" 3 (|Gastaldello et al.ll2007M . 
As can be seen (lower panel Fig [T] t ~ 0.8 — 2 Gyr), the 
group is similarly stabilized in the CC state with small en- 
tropy oscillations. 

We have done our calculations with different levels of 
AGN heating and thermal conduction, and found that our 
results are quite robust. Higher levels of AGN heating usu- 
ally correspond to larger amplitude entropy oscillations in 
the CC state. We also consider a model (run A1795-3) where 
AGN feedback only triggers conduction without heating the 
ICM (i.e., -Eagn = erg), and find that while the cluster 
cools significantly to low central entropy, it does not end up 
in a cooling catastrophe, and instead eventually ends up as 
a CC cluster as well (see the short-dashed line in Fig. [TJi 
when t < 6.4 Gyr). Note that conduction is regulated by 
AGN feedback in this run before t ~ 6.4 Gyr, after which it 
is fixed to be / = 0.4 without regulation (see further discus- 
sion in 333}. 

Although present in above runs (to account for the shut- 
off of conductivity by the HBI), the regulation of conductiv- 
ity is not required for the stability of the CC state. In run 
A1795-4, we considered a model where conductivity is trig- 
gered and then fixed to be / = 0.2 without regulation since 
t ~ 4 Gyr; the fixed lower value is meant to simulate a sit- 
uation where the field line geometry reaches a steady state 
balance between the HBI and some other mechanism (e.g, 
stirring by galaxies). The cluster evolution is similar to other 
runs. Note that for this lower value of /, the cluster would 
reach a cooling catatrophe without AGN feedback heating. 
In addition, for a fixed value of /, lower temperature clus- 
ters, and galaxy groups usually reach cooling catastrophes 
if only thermal conduction operates. 

Note that a minimum amount of thermal conduction 
is usually required in our calculations, since AGN "efferves- 
cent heating" tends to be excessively centrally concentrated: 
the very central regions of the cluster are overheated, while 
outer regions develop a cooling catastrophe. Furthermore, 
thermal conduction is the dominant energy source in mas- 
sive high-temperature clusters; without it, a much larger 
£ agn (inconsistent with observations) is required to main- 
tain the ICM in the CC state. On the other hand, in low 
temperature groups such as NGC 4325, AGN heating alone 
suffices, as seen in the lower panel of Fig. [1] 

3.2 The cycle between CC and NCC states 

Stron g AGN outbursts, e.g., H ercules A dNulsen et all 
2005h and MS0735 6+7421 (|McNamara et all 120051 
McNamara et al J 12003 ). with i? ag n up to ~ 10 62 ergs have 
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Figure 2. Time sequence of entropy in a typical model for the 
cluster Abell 1795 (run A1795-1 from t = 4 Gyr to 10 Gyr; top 
panel) and for the group NGC 4325 (run NGC4325-1 from t = 1 
Gyr to 5 Gyr; bottom panel). Each line in the top panel is plotted 
every 0.1 Gyr and that in the bottom panel is plotted every 0.05 
Gyr. The cluster A1795 oscillates in the CC state (lower branch) 
due to sporadic AGN feedback heating before t ~ 6.1 Gyr, when 
the cluster is quickly heated by a strong AGN outburst to the 
NCC state. The lines then concentrate in the NCC state (upper 
branch). The evolution of Abell 1795 is clearly bimodal, while 
bimodality in the evolution of NGC 4325 is not evident. 



been found in X-ray observations. In this subsection, we 
consider the impact of such strong AGN outbursts on the 
evolution of cool core clusters. We did our calculations for 
both the cluster A1795 and the group NGC 4325, aiming 
to investigate if such strong AGN outbursts can heat the 
CC system into a NCC state and what the role of thermal 
conduction is in this context. 

In run A 1795-1, the first AGN outburst after t = 6 
Gyr is modified to be much stronger (£7 a gn = 3 x 10 61 erg). 
Figure QJa) clearly shows that the central gas entropy is 
quickly boosted from ~ 15 keV cm 2 to ~ 130 keV cm 2 at 
t ~ 6.1 Gyr. The cluster indeed reaches a NCC state; the 
radial entropy profile is shown in the upper branch of Fig- 
ure a). During the strong AGN outburst, AGN heating 
is much stronger (by one order of magnitude) than conduc- 
tive heating. However, given the spatial dependence of AGN 
heating, conduction is still important in transporting energy 
within the cluster to offset cooling in certain regions. Since 
the cluster's temperature profile in the NCC state is nearly 
isothermal, the HBI timescale, which scales as (dlnT/dr) 1 ^ 2 
(Quatacrt 2008), is very long. We therefore assume that for 
run A1795-1, thermal conduction in the NCC state does not 
decay. Figure [2 a) clearly shows that the ICM adjusts it- 
self to the NCC state where cooling is balanced by thermal 
conduction alone. Global stability of such conduction only 
models is possible in the hottest clusters for a high level of 
conductivity (GOR08). 



Run A1795-3 is a model where AGN feedback only 
triggers conduction without heating the ICM directly (i.e., 
Sagn = erg) . When t < 6.4 Gyr, the cluster reaches and is 
then maintained in the CC state by AGN-regulated conduc- 
tion (see 33.1)1 . After t ~ 6.4 Gyr, we assume that AGN is 
continuously active so that conduction with / = 0.4 does not 
decay with time. The short-dashed line in Figure ([TJa clearly 
shows that the massive cluster A1795 is also heated to the 
NCC state, though this process takes much longer than that 
in run A1795-1. This shows that in the hottest clusters, if 
AGN can trigger a high level of conductivity, their heat in- 
put is unimportant. This is not true if the triggered level of 
conductivity is lower. In run A1795-4, we assume that con- 
ductivity is triggered and then fixed to be / = 0.2 since t ~ 4 
Gyr. The cluster evolution is similar to other runs, with the 
strong AGN outburst bringing the cluster to a NCC state. 
However, if the AGN contributes negligible heat input, as in 
run A1795-3, the lower level of conduction in this run would 
be unable to prevent a cooling catastrophe. 

Since thermal conduction will eventually decay due to 
the HBI, the cluster is expected to cool from the NCC state 
to the CC state. The HBI linear growth time for the CC clus- 
ter A1795 is t H Bi ~ 1.0 x 10 s yr; however, in the NCC state, 
the temperature profile is very flat and thus the H BI growth 
time, which scales as (dlnT/dr) 1 ^ 2 (Quatacrt 2008), becomes 
much longer. Furthermore, simulations show that the cluster 
takes several instability growth times for the magnetic field 
lines to be appreciably re-oriented. Thus, in run A1795-2, we 
assume that conduction decays on a timescale of £hbi = 1 
Gyr after the strong AGN outburst is turned off. The dot- 
ted line in Figure fTJa demonstrates that the cluster cools 
gradually until triggering new AGN activity and then stays 
in the CC state. Thus, the cluster cycles between the CC 
state and NCC state due to strong AGN outbursts and the 
HBI. 

For the group NGC 4325, we modify the AGN outburst 
at t ~ 1.9 Gyr to be much stronger (-E agn = 2 x 10 59 erg) 
and found that the group is heated to the NCC state. Fig- 
ure [Tp clearly shows that the group cycles between the CC 
state and NCC state for all three models considered in this 
paper (see Table [T] for parameters in each model), with the 
key difference that unlike the cluster case, the group does 
not stay in the NCC state for long. We have also done our 
calculations for the cluster Abell 2199 (T ou t = 4.6 keV) and 
found that a strong AGN outburst with i? ag n = 5 x 10 60 erg 
is able to heat the CC cluster to its NCC state. Thus , our 
calculations suggest the trend that more massive systems 
require stronger AGN outbursts to reach the NCC state. 

From Figure (JXJ , a scenario for the evolution of galaxy 
clusters and groups may arise naturally: a cool core group 
or cluster is maintained by normal AGN feedback heating, 
and may be heated by a strong AGN outburst to the NCC 
state, from which the system may gradually cool to the CC 
state again due to the decay of conductivity by the HBI. 
The strong AGN outburst may be triggered by a sudden 
increase in close gas supply, or mergerqj. This scenario nat- 

2 Indeed, little is know about the mechanical variability of AGN; 
possible reasonable assumptions are that the luminosity has 
a log-normal distributi on with a 'flicker-noise' power spectrum 
llNipoti fc B inncv 2003). If our hypothesis is correct, then demog- 
raphy, in particular the relative fraction of CC and NCC systems, 
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Figure 3. Histogram of gas entropy at O.Olrsoo for a sample of 28 
nearby galaxy groups from the Two- Dimensional XMM-Newton 
Group Survery jjohnson et al. 2009a). If the selection effect is not 
important, the histogram suggests that the group distribution is 
unimodal with only one peak for CC groups. 



urally explains current X-ray observations of both CC and 
NCC groups and clusters. The typical timescale of the HBI 
is around 0.1 Gyr for the cool core state of A1795 and varies 
as tuBi oc (dlnT/dr) 1/2 (|Quataertl 2008) during the cluster 
evolution. Clusters with flatter temperature profiles usually 
have larger £hbi, and thus clusters with quite flat tempera- 
ture profiles may dominate in the population of NCC clus- 
ters. Thi s agrees with current X -ray cluster observations (see 
Fig. 6 of lSanderson et alj|200fj). 

Calculations by iMcCarthv et all ((2008) suggest that 
heating a pure cooling flow cluster to a NCC cluster requires 



extremely large amounts of energy (~ 10 erg). In our cal- 
culations, AGN feedback is triggered long before strong cool- 
ing flows form, and thus the needed AGN energy is much 
less. We have performed simulations where AGN outbursts 
are triggered when a strong mass inflow forms ("the cool- 
ing flow state"), and found that much larger (several times) 
AGN heating is required to heat the cluster to its NCC state. 
The absence of strong cooling flows in clusters suggests that 
the cluster stabilizes at a CC state far before strong cool- 
ing flows form. This permits a NCC state to be attained by 
heating from a CC state, as we have seen in our simulations. 



3.3 The bimodality of CC and NCC states 

Both CC and NCC groups and clusters ha v e bee n ob- 
served in nature. Rec ently, ISanderson et al.l (|2009h and 
ICavagnolo et all (|2009h demonstrate that the CC/NCC bi- 
modality does exist in clusters: Figure 12 in the former shows 
the bimodality in the distribution of the logarithmic slopes 
of radial entropy profiles, while Figure 6 and 7 in the latter 
demonstrate the bimodality in the distribution of central gas 
entropy. 

In Figure[^a), radial entropy profiles in run A1795-1 are 
plotted every 0.1 Gyr since t = 4 Gyr, when the ICM reaches 
the CC state (see Fig. [1J,). As clearly seen, the ICM goes 
through minor oscillations in the CC state before t ~ 6.1 
Gyr (lines concentrate in the lower branch) . When the strong 



may hold the key to understanding the frequency of outbursts as 
a function of energy. 



AGN outburst is triggered at t ~ 6.1 Gyr, the cluster is then 
quickly heated to the NCC state, and stays there due to 
conductive heating in the NCC state: the lines concentrate 
there (upper branch). We clearly see the bimodality in the 
cluster evolution. Thus our simulations suggest that strong 
AGN outbursts may heat the CC cluster to its NCC state 
and the cluster CC/NCC bimodality naturally appears in 
this scenario. 

In run A1795-1, conduction is regulated by AGN feed- 
back and the HBI. In contrast, we have also considered a 
model where conduction is fixed (uncorrelated with AGN 
activity; run A1795-4). As shown in Figure [TJ a), we found 
that the cluster is also heated by the strong AGN outburst 
at t ~ 6 Gyr to the NCC state, where the cluster then stays 
due to conductive heating during our whole simulation. It 
seems that the regulation of thermal conduction is not nec- 
essary for the CC/NCC bimodality. However, note that in 
fixed-conductivity models, / must has a value within a small 
range: if / is too large, the cluster will evolve to the NCC 
state without staying in the CC state (e.g. run A1795-3 after 
t ~ 6.4 Gyr); if / is too small, the cluster may not stay in 
the NCC state, but instead develop a cooling catastrophe. 
In other words, if / is too large, we will not see a signif- 
icant population of massive CC clusters; if / is too small, 
we will not see a significant population of NCC clusters. On 
the other hand, time-varying conductivity regulated by the 
AGN could naturally circumvent this 'fine-tuning' problem: 
the cluster can stay in the CC state due to the alternation 
of radiative cooling and intermittent heating by AGN feed- 
back and conduction; the cluster can also stay in the NCC 
state, where efficient conduction triggered by strong AGN 
outbursts offsets radiative cooling. Furthermore, since con- 
ductivity increases strongly with temperature, higher tem- 
perature clusters may stay in the NCC state for a much 
longer time, which is consistent with the observational find- 
ing that the fraction of NCC systems i n clusters increases 
with the cluster mass |Chen et alj|2007l ). 

Figure EJb) shows the same plot, but for the 1 keV 
group NGC 4325 (run NGC4325-1). The lines are shown 
every 0.05 Gyr from t = 1 Gyr, when the ICM is in the 
CC state (see Fig. QJ>) to t = 5 Gyr. At t ~ 1.9 Gyr, the 
group is heated by a strong AGN outburst to the NCC state. 
However, since the conductive heating is inefficient in low- 
temperature systems (« oc T 5/2 ), the ICM can not be main- 
tained in the NCC state due to radiative cooling and cools 
to the CC state again (also see Fig. [TJ. In run NGC4325-3, 
conductivity triggered by the strong AGN outburst is mod- 
ified to be as large as the Spitzer value (see Table 1) and we 
still found that the group evolution is similar to that in run 
NGC4325-1 (see Fig.QJ>). We tried to build an equilibrium 
model with conduction alone for the NCC group NGC 4325 
and found that conduction with 10 times Spitzer value is 
required to balance radiative cooling. Obviously, in Figure 
([2])b, the CC/NCC bimodality is not evident in the group 
evolution; instead, the group distribution is unimodal: lines 
are only concentrated in the CC state, and the group distri- 
bution in NCC states is continuous rather than peaked. If 
strong AGN outbursts are not common, NCC groups may 
be rare, since they can not be maintained in NCC states by 
conductive heating. 

To test our model, we turn to group observations in 
literature. In Ras mussen fc Ponmanl (20071 1 . 14 of their 15 
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groups observed by Chandra are in the CC state, which 
may be due to a selection effect, since all of their groups 
are reasonably X-ray brig ht. Here we adopt a sa mple of 28 
nearby groups studied bv Ijohnson et a l. (2009a), which is 
the largest sample to date with high-quality XMM-Newton 
data. There exist large Chandra samples, but these only se- 
lect bright groups. We plot the distribution of the central 
gas entropy at 0.01r 5 oo in Figure ([3}, where r 5 oo is the ra- 
dius within which the mean density of the group is 500 times 
the critical density. Figure [3] clearly shows that the central 
entropy distribution is unimodal and that groups with high 
values of central entropy (S > 20 keV cm 2 ) are rare. Thus 
current group observations seem to agree with our results. 
However, note that NCC groups are very faint, and many of 
them may have not been observed yet. If more NCC groups 
are observed in future, we predict that their distribution 
is continuous, instead of peaked around one specific NCC 
state. A large group sample carefully controlled for selection 
effects would be required to do a more reliable test of our 
model. 



4 DISCUSSION 

Let us briefly summarize our findings. From a suite of ID hy- 
drodynamic simulations, we find that clusters can cycle be- 
tween CC and NCC states, driven by time-variable conduc- 
tion and/or AGN outbursts. A strong AGN outburst com- 
bined with conduction could heat a CC group or cluster to 
the NCC state. During this transition, AGN usually provides 
most of the heating energy, while conduction is important 
in transporting energy within the cluster to offset cooling in 
certain regions, given the spatial dependence of AGN heat- 
ing. The relative importance of conduction increases with 
cluster temperature, due to the strong temperature depen- 
dence of conductive flux, F <x J 15 / 2 . High temperature clus- 
ters, provided that conduction in the 'on' state is relatively 
unsuppressed / ~ 0.4 for an extended time, can even reach 
the NCC state with no energy input from the AGN. In this 
case, the AGN may simply serve as a 'switch' to regulate 
conductivity, perhaps by straightening field lines via the pro- 
duction of rising bubbles. Lower temperature clusters (or 
high temperature clusters if the maximum value of conduc- 
tion is still relatively weak) require a combination of AGN 
and conductive heating to attain the NCC state. In both 
cases, if conduction continues to operate, the cluster can 
remain stably in the NCC state. On the other hand, if con- 
duction decays via the HBI, then the cluster will cool and 
revert back to the CC state, where it remains stably with 
normal AGN feedback ( i]3.1[) until the next strong outburst 
and/or strong increase in conductivity continues the cycle. 
At the low temperature end, groups cannot be stabilized 
by any means in the NCC state, and rapidly cool to the CC 
state until the next outburst. The duty cycle, or timescale to 
cycle between CC and NCC states, shortens with declining 
temperature. 

If this hypothesis for the origin of CC/NCC cluster 
cores is correct, a number of interesting conclusions fol- 
low. Since the stabilizing effects of conduction decline with 
temperature, the NCC/CC bimodality should be a func- 
tion of temperature, being the most sharply defined for 
high temperature clusters, and vanishing in galaxy groups. 



We show there may be tentative observational evidence for 
lack of bimod al ity in core entropy in the group sample of 
Ijohnson et al.l (|2009al ) . although selection effects have to 
be carefully quantified before one can draw firm conclu- 
sions. The relative abundance of CC/NCC clusters may 
gives insights on the duty cycle on which conduction (and/or 
AGN outbursts) vary. For instance, the roughly equal frac- 
tion of CC/NCC clusters suggests that the AGN duty cy- 
cle between strong outbursts in CC systems is of order the 
HBI/cooling timescale in the corresponding NCC systems, 
though a more detailed study attempting to match the frac- 
tion of time a cluster spends at a given central entropy, to 
the distribution of entropies in the cluster population as a 
whole, would be interesting. Once again, groups can per- 
form a critical test, since they will not spend much time 
in a high-entropy state. Furthermore, since the duty-cycle 
between strong AGN outbursts is shorter in groups, tur- 
bulence or convective effects due to AGN activity which 
leave an imprint on the metallicity or entropy profile might 
have a more pronounced effect there. One possible exam- 
ple is the turbulent diffusion of metals l|Rebusco et al.ll2005l ; 
Sh arma et ail 120091) . In addition to heating the CC clus- 
ter to the NCC state, strong AGN outbursts could po- 
tentially remove the centrally-peaked metalicity distribu- 
tion observed in CC systems, resulting in a relatively flat 
metallicity profile in the NCC sta t e. W hile this distinction 



was seen m 



De Grandi fc Molendil (l200ll). more recent stud- 
ies bvlBaldi et al l <|2007|) .I Leccardi fc Molendil (2008), and 
ISanderson et alj (|2009h found that outside the very inner- 
most regions, metallicity profiles were consistent with a sin- 
gle power-law at all radii for both CC and NCC clusters. 
In contrast, the metallicity profi les in NCC groups are much 
flatter than those in CC groups l| Johnson et al.ll2009bl ). This 
is consistent with the model presented in this paper. In par- 
ticular, NCC clusters can be stabilized by thermal conduc- 
tion for sufficiently long periods for the metallicity gradient 
to be re-established, while NCC groups have a much more 
recent origin due to the short duty cycle of strong AGN 
outbursts, and thus retain evidence of AGN 'stirring' in the 
metallicity profile |f] 

Our ID calculations are frankly exploratory and simpli- 
fied in nature. Our hope is to show that AGN outbursts and 
time-variable conductivity are a plausible means of regulat- 
ing the bimodality between NCC and CC systems, motivat- 
ing future, more detailed work. Of course, the greatest gap 
in our understanding is the actual means by which the HBI 
can be counteracted to allow thermal conduction to operate. 
Even if the rising bubbles do not cause a radial reorienta- 
tion of field lines, as we have suggested, some mechanism 
(perhaps stirring of the gas by galaxies or subhalos) must 
be counteracting the onset of the HBI in clusters; otherwise, 
conduction is no longer a viable heating mechanism. This 
alone would require significant revision of theoretical mod- 
els, since no known heating mechanisms (such as AGN heat- 
ing or dynamical friction) acting alone without conduction is 
sufficient to offset a cooling catast rophe in massive clusters 
(e.g, see IConrov fc Ostrikerl [2008): such mechanisms tend 
to be too centrally concentrated toward the core, and only 
marginally sufficient energetically. Furthermore, it would 
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Table 1. List of Simulations. 



Scmt^ -^agi^ tHBl^i E agn,s fs *HBI,s 

Run keV cm 2 (10 60 erg) (Gyr) (10 60 erg) (Gyr) 



A1795-1 


15 


2.0 


0.4 





30 


0.4 


CO 


A1795-2 


15 


2.0 


0.4 





30 


0,1 


1 


A1795-3 


15 


0.0 


0.4 








0.1 


CO 


A1795-4 


15 


2.0 


0.2 


CO 


30 


0.2 


CO 


NGC4325-1 


5 


0.03 


0.4 





0.2 


0.1 


CO 


NGC4325-2 


5 


0.03 


0.4 





0.2 


0.4 


1 


NGC4325-3 


5 


0.03 


0.4 





0.2 


1.0 


CO 



a Each AGN heating episode is triggered when gas entropy (S = k^T/n c ) at r = 5 kpc drops below S CT it- 

^ The mechanical energy released during a weak (i? a gn) or strong (-E a gn,s) AGN outburst. We assume that each AGN 

outburst heats the ICM for a duration of i agn = 1-0 X 10 8 yrs. 

c The conduction suppression factor relative to the Spitzer value when AGN heating (/) or a strong AGN outburst (f s ) is 
active. 

^ Conduction is on during each active AGN heating episode. When AGN is turned off, conductivity then decays expo- 
nentially in a timescale of *hbi (after weak AGN outbursts) or tjjBl s (after strong AGN outbursts). *hbi = 00 indicates 
non-decaying conduction, while tHBl = indicates that conduction is turned off once AGN heating is shut off. 



then seem a remarkable coincidence that if one simply uses 
observed temperature profiles in clusters to construct the 
Spitzer conductive flux from the cluster outskirts, it is very 
nearly equal to that required to balance the radiative cooling 
rate as indicated by the observed X-ray surface brightness 
profile, for some reas onable fraction / ~ 0.3 o f the Spitzer 
value (e.g., Fig. 17 of lPeterson fc Fabiadl2006h . There is no 
reason in principle why such close agreement should exist, 
and seems a tantalizing hint that nature somehow 'knows' 
about Spitzer conductivity. Much work remains to be done 
before we understand if there are large secular variations to 
the apparent thermal equilibrium in clusters. 



ACKNOWLEDGMENTS 

We thank the referee, Chris Reynolds, for a helpful report. 
FG thanks Trevor Ponman and Ewan O' Sullivan for helpful 
discussions and William Mathews for a careful reading of 
the manuscript. SPO thanks MPA for hospitality. We ac- 
knowledge support by NASA grant NNG06GH95G. 



REFERENCES 

Asai N., Fukuda N., Matsumoto R., 2007, ApJ, 663, 816 
Balbus S. A., 2000, ApJ, 534, 420 

Baldi A., Ettori S., Mazzotta P., Tozzi P., Borgani S., 2007, 
ApJ, 666, 835 

Begelman M. C, 2001, in Hibbard J. E., Rupen M., van 
Gorkom J. H., eds, ASP Conf. Ser. Vol. 240, Gas and 
Galaxy Evolution. Astron. Soc. Pac, San Francisco, p. 363 

Birzan L., Rafferty D. A., McNamara B. R., Wise M. W., 
Nulsen P. E. J., 2004, ApJ, 607, 800 

Bogdanovic T., Reynolds C. S., Balbus S. A., Parrish I. J., 
2009, preprint (arXiv: 0905.4508) 

Bower R. G, McCarthy I. G, Benson A. J., 2008, MNRAS, 
390, 1399 



Burns J. O., 1990, AJ, 99, 14 

Burns J. O., Hallman E. J., Gantner B., Motl P. M., Nor- 
man M. L., 2008, ApJ, 675, 1125 

Cavagnolo K. W., Donahue M., Voit G. M., Sun M., 2008, 
ApJ, 683, L107 

Cavagnolo K. W., Donahue M., Voit G. M., Sun M., 2009, 
ApJS, 182, 12 

Chen Y., Reiprich T. H., Bohringer H., Ikebe Y., Zhang 

Y.-Y., 2007, A&A, 466, 805 
Conroy C, Ostriker J. P., 2008, ApJ, 681, 151 
De Grandi S., Molendi S., 2001, ApJ, 551, 153 
Dursi L. J., Pfrommer C, 2008, ApJ, 677, 993 
Ettori S., Fabian A. C, Allen S. W., Johnstone R. M., 

2002, MNRAS, 331, 635 
Gastaldello F., Buote D. A., Humphrey P. J., Zappacosta 

L., Bullock J. S., Brighenti F., Mathews W. G, 2007, ApJ, 

669, 158 

Guo F., Oh S. P., 2008, MNRAS, 384, 251 
Guo F., Oh S. P., Ruszkowski M., 2008, ApJ, 688, 859 
Johnson R., Ponman T. J., Finoguenov A., 2009a, MNRAS, 
395, 1287 

Johnson R., Ponman T. J., Rasmussen J., Finoguenov A., 

2009b, submitted to MNRAS 
Kaiser C. R., Binney J., 2003, MNRAS, 338, 837 
Leccardi A., Molendi S., 2008, A&A, 487, 461 
McCarthy I. G, Babul A., Bower R. G, Balogh M. L., 

2008, MNRAS, 386, 1309 
McNamara B. R., Kazemzadeh F., Rafferty D. A., Birzan 

L., Nulsen P. E. J., Kirkpatrick C. O, Wise M. W., 2009, 

ApJ, 698, 594 

McNamara B. R., Nulsen P. E. J., Wise M. W., Rafferty 
D. A., Carilli C, Sarazin C. L., Blanton E. L., 2005, Nat, 
433, 45 

Nipoti C, Binney J., 2005, MNRAS, 361, 428 

Nulsen P. E. J., McNamara B. R., Wise M. W., David L. P., 

2005, ApJ, 628, 629 
Oh S. P., Benson A. J., 2003, MNRAS, 342, 664 
O'Neill S. M., DeYoung D. S., Jones T. W., 2009, ApJ, 694, 



© 0000 RAS, MNRAS 000, 000-000 



Could AGN Outbursts Transform Cool Core Clusters? 9 



1317 

Parrish I. J., Quataert E., 2008, ApJ, 677, L9 

Parrish I. J., Quataert E., Sharma P., 2009, preprint 

(arXiv:0905.4500) 
Parrish I. J., Stone J. M., 2005, ApJ, 633, 334 
Parrish I. J., Stone J. M., 2007, ApJ, 664, 135 
Parrish I. J., Stone J. M., Lemaster N., 2008, ApJ, 688, 905 
Peterson J. R., Fabian A. C, 2006, Phys. Rep., 427, 1 
Poole G. B., Babul A., McCarthy I. G., Sanderson A. J. R., 

Fardal M. A., 2008, MNRAS, 391, 1163 
Quataert E., 2008, ApJ, 673, 758 

Rasmussen J., Ponman T. J., 2007, MNRAS, 380, 1554 
Rebusco P., Churazov E., Bohringer H., Forman W., 2005, 

MNRAS, 359, 1041 
Ruszkowski M., Begelman M. C, 2002, ApJ, 581, 223 
Ruszkowski M., Enfilin T. A., Briiggen M., Heinz S., Pfrom- 

mer C, 2007, MNRAS, 378, 662 
Sanderson A. J. R., O'Sullivan E., Ponman T. J., 2009, 

MNRAS, 395, 764 
Sanderson A. J. R., Ponman T. J., O'Sullivan E., 2006, 

MNRAS, 372, 1496 
Scannapicco E., Oh S. P., 2004, ApJ, 608, 62 
Sharma P., Chandran B. D. G., Quataert E., Parrish I. J., 

2009, ApJ, 699, 348 
Stone J. M., Norman M. L., 1992, ApJS, 80, 753 
Voit G. M., Donahue M., 2005, ApJ, 634, 955 
Zakamska N. L., Narayan R., 2003, ApJ, 582, 162 



© 0000 RAS, MNRAS 000, 000-000 



